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We have studied the ground state properties (binding energy and charge radius) using relativistic 
mean hied formalism (RMF) for Mg-isotopes in the valley of stability to drip line region. The 
RMF densities have been analyzed in context of reaction dynamics. The calculated results of 
24-40]y|-g_|_i2Q reactions at projectile energy 240 AMeV using Glauber model with the conjunction 
of densities from relativistic mean held formalism are compared with experimental data. We found 
remarkable agreement of estimated values of reaction cross sections with experimental data except 
for ^^Mg isotope. In view of this, the halo status of ^^Mg is examined through higher magnitude 
of rms radius and small value of longitudinal momentum distribution. Finally, an effort is made to 
explore the structure of ®^Mg halo candidate using Glauber few body formalism. 
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I. INTRODUCTION 

Since last few decades, the advent in radioactive ion 
beam facilities, the drip line physics and structure of drip 
line nuclei got large attention by the nuclear science com¬ 
munity. The advancement in these facilities enrich our 
knowledge toward the new exotic phenomena such as one 
and two neutron/proton halo, bubble effect in densities, 
vanishing of shell closer effect at the drip line region and 
identification of some new magic numbers for the de¬ 
formed configuration etc. The understanding of above 
mention concepts strengthen our understanding in the 
island of inversion. One of the most exotic phenomenon 
that have been exploited extensively toward drip line is 
halo status of nuclei. The halo structure originates due to 
extremely weak bound nucleons which decouple from the 
nuclear core. The interaction cross sections of such nu¬ 
clei like and a are anomalously high 

due to large magnitude of consolidated root mean square 
(rms) radii. The one neutron/proton breakup of halo 
projectiles with stable target at near/far barrier energies 
is one of the tool to investigate the structure of such 
systems Q. The measurement of longitudinal momen¬ 
tum distribution for i 6 - 18 q isotopes after one-neutron 
breakup from isotopes using the fragment sepa¬ 

rator observed narrow value of full width half maximum 
(FWHM) for (~ 44.3 ±5.9 MeV/c), which indicates 
that is a one-neutron haloQ. In addition to this, 
the measurement of nuclear reaction cross-sections for 
19’20>22(^ P, Q show that the drip-line nucleus has a 
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halo structure. The one and two neutron removal cross 
sections and momentum distribution give indication for 
the halo status of a nuclear system. It is relevant to men¬ 
tion that 92(7 has N = 16 which is a new magic number 
in neutron-rich region i0 and forms a Borromean halo 
structure (^^C is unstable). It is of large interest to study 
the existence of bound states of two nucleons (neutrons 
or protons) with core in the Borromean structure, which 
have not appeared in vacuum. Soon after was in¬ 

cluded in the family of neutron halo. The measurement of 
interaction cross sections for Ne isotopes from the stabil- 
iW line to neutron drip line at 240 MeV/nucleon energy 
0, show the dependence of mass number for 97-32]\jg 
isotopes, which have been explained by nuclear defor¬ 
mations. The enhancements of interaction cross sections 
particularly for ^^Ne and ^^Ne have been quoted by an s- 
dominant halo structure of ^^Ne and s- or p-orbital halo 
in ^^Ne Q. The isotope "^^Ne having neutron number 
N=21, seem to break the shell closer structure. As a 
consequence, a large value of deformation associated with 
the strong intruder configuration put in at ” island of in¬ 
version” . 

Recently, M. Takechi et al. 0 , have measured the data 
of reaction cross section for Mg isotopes at energy 240 
A MeV at RIBF and RIKEN. In this study, they have 
observed large reaction cross section value for ^^Mg as 
compared to other isotopic chain and consequently pre¬ 
dicted ^^Mg to have a halo structure. In another study 
S. Watanabe concluded from their fully microscopic 
double folding frame work with the AMD densities, that 
^^Mg exhibits deformed halo structure. In our earlier 
work we explored the structure of some isotopes of Ne, 
Mg and Si [Q. The identification of this new candidate 
as halo encourage us to review our study for Mg isotopes. 
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Therefore we have made a motivated effort to study the 
reaction dynamics of Mg isotopes at energy 240 AMeV. 
A sincere effort is made to analyze the structural fea¬ 
tures of ^"^Mg. The reaction/interaction cross section, 
angular elastic differential cross section, one nucleon re¬ 
moval cross section and angular momentum distributions 
parameters are exercised for the analysis of such nuclear 
systems. 

The section II contains a brief description of Glauber for¬ 
malism. Section III describe the calculations and results, 
and finally a summary and conclusions is described in the 
section IV. 


II. THE FORMALISM 

We use the well known Glauber approach to investi¬ 
gate reaction dynamics [l3l - [T5l| . The study of reaction 
dynamics in the framework of this approach strongly de¬ 
pends on the densities of the projectile and target nuclei, 
and we have used the microscopic relativistic mean field 
(RMF) densities with NL3 parameter (T^ . 


The profile function Tnn for optical limit approximation 
is defined as 

Tnn = Tij{beff) = crNNexp{- ), (5) 

^^Pnn ^Pnn 

for finite range and 

^NN = ^ O-NN^ib), ( 6 ) 

for zero range with beff=\ b — -I- b is the im¬ 

pact parameter. Where and t are the dummy vari¬ 
ables for integration over the z-integrated target and pro¬ 
jectile densities. The parameters crNN,o:NN and Pnn 
usually depend upon the proton-proton, neutron-neutron 
and proton-neutron interactions. Here cfnn is the total 
nuclear reaction cross section of NN collision, unn is 
the ratio of the real to the imaginary part of the for¬ 
ward nucleon-nucleon scattering amplitude and Pnn is 
the slope parameter. The slope parameter determines 
the fall of the angular distribution of the N-N elastic scat¬ 
tering. 


A. Glauber Model 


1. Reaction cross section 


2. Angular elastic differential cross section 


The theoretical formalism to study the reaction cross 
sections using the Glauber approach has been given by 
R. J. Glauber [I^. The standard Glauber form for total 
reaction cross sections is expressed as [3 [11 

pOO 

aH = 2TT b[l - T{b)]db, (1) 

Jo 

where ’T(b)’ is the Transparency function with impact 
parameter ’b’. The function T(b) is calculated by 

T{b) = exp[-^cry /^( 5 )^( 1 "^ - ^|)ds]. (2) 

* j' 


Here, the summation indices i, j run over proton and neu¬ 
tron and subscript ’p’ and’t’ refers to projectile and tar¬ 
get respectively, dy is the experimental nucleon-nucleon 
reaction cross-section which depends on the energy. The 
z- integrated densities are defined as 

poo 

p{uj) = / p{\/w'^ + z'^)dz, (3) 


with Initially Glauber model was de¬ 

signed for the high energy approximation. However, it 
was found to work reasonably well for both the nucleus- 
nucleus reaction and the differential elastic cross-sections 
over a broad energy range Em. The modified trans¬ 
parency function T(b) is given by 


T{b) = exp[- f f 

JpJt ^ j 




( 4 ) 


The nucleus-nucleus elastic scattering amplitude is 
written as 


/ K f 

F{q) = — / dbe^<i \l - (7) 

27r J 

At low energy, this model is modified in order to take care 
of finite range effects in the profile function and Goulomb 
modified trajectories. The elastic scattering amplitude 
including the Goulomb interaction is expressed as 

F{q) = e^>^‘{F,oui(.q)+'-^ J 

( 8 ) 

with the Goulomb elastic scattering amplitude 

FcouM) = - 5 —exp{- 2 i? 7 ln(;^) -|- 2LargT{l + tr])}, 

IK 

( 9 ) 

where K is the momentum of projectile and q is the 
momentum transferred from the projectile to the tar¬ 
get. Here rj = ZpZ^^ jhv is the Sommerfeld param¬ 
eter, V is the incident velocity of the projectile, and 
Xs = —2r]ln{2Ka) with a being a screening radius. The 
elastic differential cross section is given by 


da 


da 

dap 


- \m\^- 

( 10 ) 

\Fiqr 

\Fcoui{qW' 

( 11 ) 
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3. one nucleon removal cross section 


B. Relativistic mean field formalism 


The one nucleon removal cross section, a-N, may be 
defined as 


U-N = Sc / dka, 


a—{k,g—0),c 


( 12 ) 


we assume that the core remains in its ground state, g=0. 
The one nucleon removal cross section consists of both 
elastic and inelastic part and can be calculated by 


^ ^el I ^inel 

O^-N — “T ^-N • 


(13) 


The cross section due to the elastic breakup process is 
given by 

While the cross section from inelastic breakup is given 

by 




_^-2ImxcT(bc)-2ImXNTibc+s) 


l</'o)|}. (15) 


4- Longitudinal momentum distribution 

The momentum distribution of core after the inelastic 
breakup of Projectile reads as: 


da^44-N 

dP 


The relativistic mean field approach is well docu¬ 
mented in Refs. [Tsl - I^ . The basic ingredient of RMF 
model is the relativistic Lagrangian density for a nucleon- 
meson many body system which is defined as 

I 22 I 3 I 4 T I 

~ 2 '^cr^ ~ 3520 - - - gs'tpi'^iCr 








(19) 


[ + fiK)\F(K,o)A<l)\X^^) 

c#0 ^ 


The Scattering wave function of nucleon is approximated 
by a plane wave reduced to 


dFii 




1 1 
{2TrH)^ 2j + 1 




,( 17 ) 


Here a, and are the fields for cr-, w- and p-meson 
respectively. is the electromagnetic field. The rpi are 
the Dirac spinors for the nucleons whose third compo- 

2 -1 

nent of isospin is denoted by . gsi gu>i gp and ^ 
are the coupling constants for the linear term of cr-, w- 
and p-mesons and photons respectively. p 2 and ps are 
the parameters for the non-linear terms of the cr-meson. 
M, rria, rriuj and rrip are the masses of the nucleons, cr-, 
w- and p-mesons, respectively, and are 

the field tensors for the V^, and the photon fields, 
respectively. The quadrupole moment deformation pa¬ 
rameter /32, root mean square radii and binding energy 
are evaluated using the standard relations [l^. The nu¬ 
clear density p = is obtained by solving the 

equation of motion obtained from the above Lagran gian 
The values of the parameters for NL3 are given as [1^ 
Ps=10.217, 5,^=12.868, 5p=4.574, 52 = -10.431 (/m"!), 
53 = -28.885, and M=939, m,^=508.194, m,^=782.501, 
TOp=763.0 in MeV. 


III. CALCULATIONS AND DISCUSSION 


where bjxr stands for the impact parameter of valence nu¬ 
cleon with respect to the target. The longitudinal mo¬ 
mentum distribution obtained by the integration of above 
equation over transverse component of momentum Fj_, 
gives. 





da^d4-M 

dP 


1 

2iTh 


X 


J db^il - J J^g^-2ImXCT{bM-s) 

dze^F\^Unij{r)Yimi{f)\'^. (18) 


(2Z + 1) 


The ground state properties of Mg-isotopes are calcu¬ 
lated using microscopic relativistic mean field formalism 
Fig 1 shows the calculated values of binding energy (B. 
E.) and charge radius (rc) of considered set of isotopes of 
Mg. The lower panel of the figure shows the calculated 
values of B.E, which find nice comparison with experi¬ 
mental data which looks nice agreement with each oth¬ 
ers. A deep inspection of the figure shows that the B. E.’s 
slightly under estimate up to the ^®Mg and overestimate 
beyond ^^Mg isotopes. While the upper panel of figure 
compare the calculated value of rms charge radius {tc) 
of considered set of isotopes with the experimental data. 
The calculated values of are underestimated in com¬ 
parison to experimental data. The success of relativistic 









4 



24 26 28 30 32 34 36 38 40 

Mass Number (A) 


FIG. 1: The values of B.E in MeV and Charge radius of 
Mg isotopes after obtained from RMF(NL3) as a function of 
mass number (A). The experimental data are also given for 
comparison whichever be available. 

mean field formalism depends on the appropriate choice 
of densities and we have used RMF(NL3) densities as an 
input of Glauber model. 

In the measurement of reaction parameter through 
Glauber formalism, one of the input for evaluation of 
profile function in Glauber model is its energy as well 
as isospin dependent parameters. The values of these 
parameter at Eproj= 240 AMeV are aNN= 3.266868 
(/m2), aArAr=0.6800303 and Pnn= 0.097843707 (/m^). 
These values have been estimated by spline interpolation 
from Ref. [1^ . Other important inputs of Glauber model 
are the densities of the projectile and the target nuclei. 
Fig. 2 shows the relativistic mean filed (RMF) densities 
of projectiles with NL3 parameter set, which shows the 
density profile as a function of radial distance for 
nuclei. Here the nucleonic densities distribution are of 
larger values at the center and goes on decreasing as the 
radius increases, where the small depletion in densities 
also appear at the center for these nuclei. One may also 
observed from the figure that the skin effect increases 
with increase of isotopic mass number. 

These densities can be feed as an input of Glauber model 
after converting into spherical equivalent of it in terms of 
Gaussian coefficients. We have converted these densities 
into Gaussian form and calculated their values in terms 
of Gaussian coefficients Ci’s and a^’s using relation: 

n 

= X! c^exp[-a^r'^], (20) 

i=l 

The Gaussian coefficients which are used as input in the 
Glauber model code M are listed in Table 1 for NL3 
interaction. 

For the address of reaction dynamics, the single par¬ 
ticle wave function is used in Glauber model. The ra¬ 
dial part of single particle wave function have been ob¬ 
tained after solving Schrodinger equation using Wood- 
Saxon type potential as in the form: 

U{r) = -vofir) + Vis{l.s)rl^^ + VcouU (21) 



FIG. 2: The density profile of Mg isotopes from RMF(NL3) 
formalism as a function of radial distance. 


TABLE I: The Gaussian coefficients for the Projectile and 
Target after fitting RMF(NL3) densities are listed in Table 
below. 




RMF(NL3) 


Nuclei 

Cl 

fli 

C2 

02 


-0.232654 

0.638687 

0.517232 

0.339911 


-3.6936 

0.286163 

3.80452 

0.261585 


-3.89231 

0.281947 

4.00019 

0.258213 

^^Mg 

-4.07533 

0.277958 

4.18027 

0.25493 

^^Mg 

-4.065 

0.266396 

4.1637 

0.244635 

^^Mg 

-3.93233 

0.256062 

4.02441 

0.234787 

^^Mg 

-3.89795 

0.246598 

3.98355 

0.226164 

^°Mg 

-3.87287 

0.237936 

3.95204 

0.218274 

^^Mg 

-3.41668 

0.224842 

3.52235 

0.206223 

^^Mg 

-2.96194 

0.211646 

3.09244 

0.194095 

^^Mg 

-2.96276 

0.206741 

3.09102 

0.189597 

^^Mg 

-2.96402 

0.202127 

3.08998 

0.185357 

^^Mg 

-2.96363 

0.197607 

3.08718 

0.18121 

^^Mg 

-2.96752 

0.193359 

3.08865 

0.177322 

^^Mg 

-2.96752 

0.189208 

3.08637 

0.173515 

CO 

00 

-2.97189 

0.185226 

3.08797 

0.169881 

^^Mg 

-2.97547 

0.181501 

3.08903 

0.166466 

'^°Mg 

-2.98107 

0.17785 

3.09205 

0.163137 


where /(r) = [1 -I- exp(^^^)]“^ and R = The 

first term of equation (21) contains the central potential, 
second term contains spin orbital part and the last term 
of the equation contains Coulomb part of potential. A 
be the mass number of nucleus. We fixed the value of 
To = 1.2 fm and diffuseness parameter ”a” as 0.6 fm in 
our calculations. 

Fig. 3 represents the values of ap for reac¬ 

tions at iJproj =240 AMeV as a function of A of projectile 
nucleus. The calculated values of ap in the figure show a 
remarkable agreement with the experimental data except 
for the case of ^^Mg, which support the success of RMF 
densities for the study of reaction dynamics. One may see 
some deviation particularly for higher Mg isotopes and 
^^Mg projectile. This difference in theoretical and exper¬ 
imental data for ^^Mg nucleon is of further significance 
which seem to suggest that, it exhibits unusual structure 
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Mass number (A) 

FIG. 3: Reaction cross section for as projectile with 

target nucleus at Eproj= 240 AMeV. The experimental 
data are also given for comparison 0. 



0^^(deg.) 

FIG. 4: The values of angular elastic differential cross section 
for reactions at Eproj= 240 AMeV. 

and hence needs further investigation. This point is ex¬ 
plored in further discussion. 

Fig. 4 shows the angular elastic differential cross 
sections of projectiles on the carbon target at 

Eproj=2^0 AMeV. The inspection of the figure suggest 
that, two dip positions are observed at an angles 0c.m 
= 2° and 8°. One may also observe that the large dip 
appeared for ^^Mg and ^®Mg at these angles. Specifically 
speaking, the largest dip at ^^Mg projectile may be 
associated with the loosely bound structure of this 
nucleus. The above observations give indication of halo 
behavior of ^^Mg. 

For further study of the case of ^^Mg projec¬ 
tile, we have used Glauber formalism with two body 
(cored-nucleon) system. Further details of the calcula¬ 
tions can be seen in [2^, [2^. Fig. 5 represents the cal¬ 
culated values of ap (in mb) and (in fm) for the 

^^Mg projectile. The upper panel of Fig. 5 presents the 
root mean square radius r^ms and lower panel represents 
(Jr for ^"^Mg as a projectile over the target at energy 
240 AMeV as a function of diffuseness parameter ’a’ in 
fm. We have tried to fit the reaction cross section for 
^^Mg projectile at different values of diffuseness param¬ 



a (fm) 

FIG. 5: (A) Upper panel of the figure show the rms radius 
(frms) in fm and (B) lower panel show the values of reaction 
cross section (jr in mb as a function diffuseness parameter ’a’ 
in fm for ®^Mg. 



r (fm) r (fm) 

FIG. 6: (A) Comparison of single particle wave function for 
different values of diffuseness parameter ’a’ and (B) Gompar- 
ison of the Wood Saxon densities with different diffusion pa¬ 
rameter and RMF(NL3) density of ^^Mg as function of radial 
distance in fm. 


eter. It is clear from the figure that the reaction cross 
section and diffuseness parameter are linearly dependant 
to each other. The value of ap obtained at diffuseness 
parameter = 3.5 fm is 1537 mb, which is well comparable 
to the experimental observation 1536±15 mb. Hence the 
lower panel of the figure shows that the value of reaction 
cross section fit with the experimental value at a = 3.5 
fm. The upper panel of the figure shows the rms radius 
values of core-l-neutron system as a function of diffuse¬ 
ness parameter. We observed that the value of root mean 
square radius of projectile (core -I- nucleon) at a=3.5 fm 
is 6.041 fm. Thus the large value of reaction cross sec¬ 
tion has direct consequence to their large radius and halo 
nature of ^’^Mg isotopes. 

Fig.6 represents the comparisons of wave function for 
different value of diffuseness parameters and RMF den¬ 
sity for ^’^Mg with Wood-Saxon density having diffuse¬ 
ness parameter 0.6 fm and 3.5 fm. The left side of the 
figure show variation of single particle wave function with 
radial distance for different values of diffuseness param¬ 
eter. This figure suggest that the wave function is more 
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FIG. 7: Variation of Total reaction cross sections for as 
function of Eproj- 



Ep.,(AMeV) 

FIG. 8: One neutron removal cross section for the reaction 
reaction including elastic and inelastic part as a 
function of Eproj A MeV 


steeper at diffuseness value of 0.6 fm where as more 
broader at 3.5 fm. The right side of this figure show 
that the behavior of Wood-Saxon density is similar to 
RMF at a=0.6 fm. At the value of ’a = 3.5 fm’, we 
get a long extension of density even beyond ^15 fm, 
further suggesting the possibility of a halo structure of 



FIG. 9: Longitudinal momentum distribution of ®®Mg from 
the ^^Mg+^’^C at projectile energy 240 A MeV. 


^"^Mg nucleus. Although, this large value of a = 3.5 fm 
looks unusual, but still be relevant for exotic nuclei ex¬ 
hibiting halo structure. For further significance of higher 
diffuseness parameter ” a”, one need to have extensive 
investigation by addressing reaction dynamics involving 
various halo nuclei. We have fitted density with a Wood- 
Saxon form taking variation in diffuseness parameter ’a’. 
Then, we used these densities in the reaction calculation 
to evaluate the reaction cross-section <jr. We find that 
the fitted density reproduce the experimental data for 
^^Mg at a = 3.5 fm. 

Fig. 7 represents the variation of reaction cross section as 
a function of projectile energy for 30-1000 AMeV using 
wave function with a=3.5. The experimental values are 
also given for comparison. Fig. 8 shows the one neutron 
removal cross section for ^^Mg-|-^^C reaction as a function 
of Eproj- One neutron removal cross section consists of 
both elastic and inelastic component. It is clear from the 
figure that the inelastic component in neutron removal 
cross section dominate over their elastic component at 
Eproj=24:0 A MeV. The trend of reaction cross section 
in figure 7 and one neutron removal cross section in figure 
8 are similar, but small hike is appeared in one neutron 
removal cross section at Eproj=l00 A MeV, because of 
its elastic component. 

Fig. 9 show the calculated longitudinal momentum dis¬ 
tribution of ^®Mg core for the reaction ^^Mg-|-^^C at 
Eproj=2A0 A MeV. The trend of distribution exhibits 
the Gaussian pattern. So we compare our calculated val¬ 
ues of longitudinal momentum distribution using RMF 
densities with in Gaussian function. The circle points 
show the calculated values of momentum distribution of 
one nucleon from the ^^Mg projectile and solid line is the 
fitted gaussian curve. By fitting curve the observed value 
of FWHM comes out to be 80.02 MeV/c. 


IV. SUMMARY 


In summary, we have calculated the ground state prop¬ 
erties of Mg isotopes and also studied the reaction cross 
sections of these isotopes taken as projectile from the 
valley of stability to drip line region with stable tar¬ 
get at Eproj 240 AMeV. We found remarkable agreement 
of ground state properties of Mg-isotopes with available 
data. The estimated values of reaction cross section using 
densities from RMF formalism are nicely compered with 
the experimental data. The excellent agreement of esti¬ 
mated reaction cross section values except for ^"^Mg iso¬ 
tope is an evidence of predictive power of RMF. The skin 
effect variation with mass number is studied in context of 
density profile. The study of angular elastic differential 
cross section for 34-38]y[g further investigation with 
Glauber two body calculation also support the halo sta¬ 
tus of ^"^Mg. Subsequently we examined the halo status 
of ^^Mg and it seems justified from its higher magnitude 
of rms radius ~6.041 fm and small value of FWHM (80.02 
MeV/c) of longitudinal momentum distribution. 
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